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CHAPTER  I 


INTRODUCTION 


1 . 1  Background 

When  an  incident  sound  field  encounters  an  obstacle,  the 
acoustic  energy  is  scattered  in  all  directions.  The  scattered 
acoustic  field  depends  on  geometrical  considerations  such  as 
curvature,  size,  discontinuities  in  slope  (sharp  edges)  and  in 
curvature  as  well  as  the  impedance  of  the  surface  or  its  elasticity. 
If  the  incident  acoustic  field  is  harmonic  in  time,  then  the 
scattering  field  also  depends  on  the  wavelength  vs.  a  typical  length 
dimension  and  the  scattered  field  generates  a  directivity  that  is 
fixed  in  time  and  is  the  result  of  interference  between  different 
parts  of  the  scatterer.  If  the  incident  field  is  an  acoustic  pulse, 
then  the  scattered  field  is  time  dependent  and  the  directivity  is 
time  dependent  and  is  the  result  of  distinct  time  signatures  of  each 
part  of  the  scatterer.  If  the  surface  impedance  is  frequency 
dependent,  then  the  scattered  field  is  influenced  by  the  interference 
between  the  frequency  of  the  incident  field  and  that  of  the  impedance 
of  the  surface  of  the  scatterer.  However,  the  impedance  of  the 
scatterer  would  not  have  the  same  interference  with  an  incident 
acoustic  pulse.  This  difference  between  continuous  and  pulse 
acoustic  scattering  motivated  this  study. 

1 . 2  Basic  Approach 

To  investigate  the  influence  of  surface  impedance  on  the 
reflection  of  an  acoustic  pulsed  wave,  a  plane  pulsed  wave  reflecting 
from  an  impedance  covered  plane  is  investigated.  This  geometry's 
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simplicity  allows  an  extensive  investigation  of  the  pulse  shape  and 
time  signature  on  the  reflected  pulse.  For  such  a  flat  surface,  the 
reflected  field  is  primarily  due  to  geometric-optics.  To  investigate 
the  influence  of  curvature,  the  investigation  is  extended  to  include 
the  pulse  reflection  from  an  impedance  covered  cylinder.  Integral 
transforms  on  time  and  then  use  of  eigenfunction  expansions  result 
in  a  series.  This  leads  to  the  use  of  the  Sommerf eld-Watson 
transformation  which  results  in  a  double  integral  for  the  time 
signature  of  the  reflected  field.  The  resulting  integral  is  then 
separated  into  two  parts,  i.e.,  the  geometrically  diffracted  field 
and  the  creeping  wave  field.  These  are  separately  evaluated  by 
asymptotic  methods  and  by  use  of  residue  methods  resulting  in  fastly 
convergent  series.  Those  series  are  then  integrated  in  a 
Duhamel-type  integral  to  obtain  the  scattered  field  due  to  any  type 
of  incident  pulse  time  signature. 

1.3  Review  of  Previous  Investigations 

Studies  of  the  scattering  from  objects  have  an  extensive 
history  in  acoustics  and  optics.  By  far  the  greatest  attention  has 
been  directed  to  objects  with  either  rigid  or  pressure  release 
surfaces.  Acoustic  scattering  studies  have  mainly  focused  on 
scatterers  which  are  simple,  e.g.,  spherical,  circular  cylindrical, 
elliptical,  etc.,  and  for  which  the  method  of  separation  of  variables 
can  be  brought  to  bear.  In  addition,  there  is  a  large  literature  of 
alternative  methods  that  range  from  sophisticated  function  theoretic 
techniques  to  purely  numerical  ones  in  which  integral  equations  for 
unknown  surface  distributions  are  approximated  by  sets  of  linear 
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equations  and  solved  numerically.  Keller  and  Levy1  have  analyzed  the 
diffraction  of  time  harmonic  waves  by  smooth  objects.  Using  the 
geometrical  theory  of  diffraction,  based  on  the  postulate  that  fields 
propagate  along  rays,  they  solve  the  problen  of  scattering  for  simple 
shapes  like  spheres  and  cylinders.  Keller  formulates  his  solution 
for  different  types  of  boundary  conditions,  where  he  expands  the 
solution  asymptotically  for  large  Ka.  An  important  advantage  of  this 
method  is  that  it  does  not  depend  on  separation  of  variables  or  any 
such  similar  procedure.  Therefore  the  shapes  of  the  objects  to  which 
it  can  be  applied  are  not  restricted,  but  it  can  be  expected  to  yield 
good  results  only  when  the  wavelength  is  small  compared  to  typical 
obstacle's  dimensions  and  therefore  cannot  be  applied  in  the  case  of 
a  pulse. 

Blank  and  Keller2  discuss  the  diffraction  and  reflection  of 
plane  pulses  by  wedges  and  corners  for  the  special  case  of  rigid  or 
pressure  release  boundary  conditions.  By  following  the  propagation 
of  the  plane  discontinuity  in  xyt-space,  they  obtain  an  explicit 
closed  form  expression  in  terms  of  elementary  functions.  But  their 
technique,  known  as  the  method  of  conical  flow,  cannot  be  applied  to 
the  more  general  case  of  impedance  boundary  condition.  Friedman  and 
Shaw3  have  studied  the  diffraction  of  pulses  by  cylindrical  obstacles 
of  arbitrary  cross-section.  They  discuss  the  two-dimensional  problem 
of  the  diffraction  of  a  plane  acoustic  shock  wave  in  terms  of  a 
characteristic  boundary  value  problem  in  four-dimensional  (x,y,z,t) 
to  which  Green's  identity  may  be  applied,  instead  of  an  initial 
boundary  value  problem.  An  integral  equation  for  the  surface  values 
of  the  pressure  is  formulated  and  approximated  by  a  set  of 


successive,  non-simultaneous  algebraic  equations  which  can  be 
solved  for  a  given  geometry.  For  the  case  of  a  square  box  with  rigid 
boundary  conditions ,  the  surface  values  of  pressure  are  solved  in  this 
manner  for  a  period  of  one  transition  time. 

Hickling  and  Means4 » 5  have  investigated  the  scattering  of  FM 
pulses  by  spherical  elastic  shells  in  water.  They  discuss  the 
problem  in  two  areas  of  practical  interest,  (a)  linear 
frequency-modulated  (FM)  pulses,  and  (b)  higher  frequencies  where  the 
wavelength  is  much  smaller  than  the  exterior  radius  of  the  shell. 
Adopting  a  statistical  approach  based  on  signal  processing,  they 
compute  a  cross-correlation  function  of  each  echo  with  its  incident 
pulse.  They  show  that  there  is  a  time  delay  between  the  first, 
second,  third,  and  later  peaks  in  the  cross-correlation  function  and 
reach  the  conclusion  that  the  first  pulse  in  the  echo  corresponds  to 
a  reflection  from  outer  surface  of  the  shell  closest  to  the  source 
and  the  remaining  pulses  are  due  primarily  to  the  elastic  response  of 
the  sphere. 

Shaw6  studied  the  diffraction  of  arbitrary  plane  acoustic  pulses 
by  obstacles  whose  boundary  offers  inertial  but  no  elastic  resistance. 
The  boundary  condition  is  expressed  in  terms  of  the  familiar  Robin-type 
(<})  +  K  3$/3n  =  0),  where  <p  is  tne  velocity  potential.  A  large  value  of 
K  Implies  a  heavy  scattering  surface  (rigid)  and  a  small  value 
corresponds  to  an  almost  free  surface.  The  method  of  solution  is  very 
much  the  same  as  that  used  in  Ref.  3.  Through  some  numerical 
approximation  he  solves  the  integral  equation  for  a  square  box  and 
finds  the  surface  values  of  pressure  for  large  K  and  surface  values  of 


velocity  for  small  K 
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Uberall7  investigates  the  multiple  echoes  that  are  returned  in  the 
scattering  of  a  single  pulse  and  analyzes  the  presence  of  creeping 
waves  or  circumferential  waves  in  the  diffraction  process.  For  the 
special  case  of  an  infinitely  long,  acoustically  soft  circular 
cylinder,  he  describes  the  creeping  waves  using  Sommerf eld-Watson 
transformation  method  of  Franz8  and  assumes  that  the  initial  plane 
sound  wave  possesses  an  arbitrary  pulse  shape.  The  expressions  for  the 
creeping  pulses  obtained  on  the  surface  of  the  scatterer  includes,  but 
go  beyond  previous  results  of  Friedlander8 ,  found  by  a  different 
method. 

••  •• 

Doolittle,  Uberall  and  Uguncius10  discuss  the  problem  of 
scattering  of  a  simple  harmonic  plane  wave  by  an  infinite  elastic 
cylinder.  They  apply  the  Sommerfeld-Watson  transformation  on  the  normal 
mode  series.  The  complex  velocities  of  the  ensuing  circumferential 
waves  are  found  by  numerical  search  for  the  zeros  of  a  3x3  determinant 
in  the  complex  plane.  Two  distinct  types  of  zeros  are  found:  (a) 
"Franz-type"  zeros,  similar  to  those  appearing  in  the  scattering  from 
rigid  cylinders.  These  waves,  propagating  with  speeds  below  the  sound 
speed  in  the  surrounding  medium,  are  heavily  attenuated  and  penetrate 
and  exit  the  cylinder's  surface  tangentially,  and  (b)  "Rayleigh-type" 
zeros,  also  found  earlier  by  Goodman  and  Grace1,1  and  Franz  and  Beckmann 
for  cylinders  of  finite  conductivity12.  The  waves  are  weakly  attenuated 

and  have  speeds  approaching  elastic  velocities. 

••  «• 

Ugincius  and  Uberall13  investigate  the  steady-state  scattering  of 
a  plane  harmonic  wave  from  an  infinite  elastic  circular  cylindrical 
shell,  as  treated  by  the  creeping  wave  method.  This  work  is  basically 
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an  extension  of  their  previous  study  of  sound  scattering  from  elastic 
cylinders,  but  new  interesting  results  were  obtained  in  the  limit  of  a 
thin  shell.  The  resulting  residue  series  is  generated  by  poles  that 
are  the  complex  zeros  of  a  6x6  determinant.  These  zeros  are  found 
numerically  by  an  extension  of  the  Newton-Raphson  method  for  complex 
functions.  Finally  it  was  found  that  besides  the  infinity  of  the  well 
known  rigid  zeros,  there  exists  a  set  of  additional  zeros  that  give 
rise  to  generalized  Rayleigh  and  Stoneley  waves.  The  behavior  of  these 
zeros  has  been  studied  as  a  function  of  the  wave  number  and  also  of  the 
shell  thickness. 

In  most  cases,  the  scattering  problem  was  solved  for  objects  with 
constant  impedance  or  elastic  boundary  conditions.  A  combination  of 
non-harmonic  incident  wave  and  frequency  dependent  impedance  covering 
will  create  new  complications  in  the  numerical  computation  of  the 


scattering  problem  which  will  be  discussed  later. 


CHAPTER  II 


REFLECTION  OF  A  PLANE  WAVE  PULSE  FROM  AN  IMPEDANCE  COVERED  PLANE 

2.1  Statement  of  the  Problem 

Consider  an  infinite  flat  surface  in  the  x-y  plane  which  bounds 
a  semi-infinite  acoustic  medium  z  >  0.  In  this  chapter  a  solution  is 
sought  for  the  reflected  acoustic  field  due  to  a  plane  pulse  incident 
on  an  infinite  plane  boundary.  The  boundary  is  covered  with  a 
material  possessing  a  locally  reacting  impedance,  which  is  frequency 
dependent  and  may  have  one  or  multiresonant  character. 

An  incident  plane  pulse  is  propagating  in  (-z)  direction,  normal 
to  the  boundary.  If  t  =  0  correspond  to  the  time  of  arrival  of  the 
pulse  at  the  boundary,  then  the  incident  wave  can  be  expressed  as: 

pinc(c)  = 

Applying  the  Fourier  transform  to  pinc(t)  one  obtains 

+00  +0O 

4> (co)  =  /  Pinc^  e  iUt  dt  =  /  G(t)  e  iuC  dt  *  1  (1) 

—GO  — OO 

Application  of  Fourier  inverse  transform  on  <j>(w)  would  lead  to 
reconstruction  of  the  incident  v.’ave 

5(t)  =  ~  /  4.(«)  eiUit  du  =—■  f  eiL)t  dt  (2) 

— OO  — oo 

The  integrand  in  Eq.  (2),  corresponding  to  a  definite  value  of  to, 
f(to)  =  ^(w)  e*-uC,  represents  a  plane  harmonic  wave.  Multiplying  f(to) 
by  the  reflection  coefficient  of  the  boundary  one  obtains  the 
contribution  to  the  reflected  field  arising  from  an  individual 
harmonic  wave.  Finally  applying  the  superposition  principle,  the 
reflected  field  due  to  an  incident  pulse  can  be  written  as  follows: 
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“T30 

h(t)  =  -s—  /  R(w)  e*Wt  dm 


The  reflection  coefficient  R(u>),  for  normal  incidence  is  expressed 


=  1  - 


where  Z(w)  represents  the  impedance  of  the  surface,  and  p  and  c 
represents  respectively  the  mass  density  and  sound  speed  of  the 
medium.  Considering  a  simple  oscillator  excited  by  a  time  harmonic 
force  F  =  f0e^U)t,  the  impedance  can  be  expressed  as  the  ratio  of 
driving  force  to  velocity 


Z(w)  =  T  =  c0  +  imu)0(—  -  ^2-)  , 
X  U 


where  oj0,  k  and  m  represents  the  natural  frequency,  the  stiffness  and 
the  mass  of  the  oscillator.  Dividing  Z(w)  by  pc  and  introducing  the 
dimensionless  constants  aj  and  bj  we  have  a  model  for  the  impedance 
of  the  surface: 


where 


Z(u>  .  ^ 

-  =  bi  +  iai  ( — - ) 

PC  1  1  U0  !d 


bi  =  —  ,  and  a  i  =  — 

1  pc  1  pc 


Finally  substituting  Eqs.  (4)  and  (5)  in  (3)  one  obtains: 
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h(t)  -  -5-  / 


iwt  ,  1  r 

e  du  -  ■=-  J 


(bi+l)+iai(- 


“o  t°° 

h(t)  -  6(c)  J 


neia)°tn 

(bj+Dn+iaj/n2-!) 


The  integral  in  Eq.  (6)  can  be  computed  by  using  the  Residue  Theorem 
and  closing  the  contour  of  integration  along  a  semi-circle  in  upper 
half  ri-plane.  The  poles  of  the  integrand  are  the  roots  of 


(b]+l)n  +  iai<n2-l)  =  0* 


It  could  be  shown  that  both  roots  of  Eq.  (7)  fall  in  the  upper  half 
of  n-plane.  For  (bj+l)  <  2a,  they  would  be  symmetric  with  respect  to 
imaginary  axis  and  for  (bj+1)  >  2aj,  both  roots  would  be  purely 
imaginary.  For  the  special  case  of  (bj+1)  =  2ai,  the  roots  coincide 
and  hence  the  pole  would  be  of  second  order.  Let  rin(n=1.2)  be  the 
roots  of  Eq.  (7)  and  qn  be  the  corresponding  residues 


qn  =  -  n  =»  1,2 

(bj+1)  +  2ia1nn 


Finally  Eq.  (6)  can  be  written  as: 


h(t)  -  6(t)  -  —  ( 2tt i )  2  qn 
"  n=l , 2 


It  should  be  observed  that  the  right  hand  side  of  Eq.  (9)  is  thus 
necessarily  a  real  quantity  as  is  expected.  The  pulse  response 
expressed  by  Eq.  (9)  could  be  substituted  in  Duhamel's  integral  to 


obtain  the  reflected  pressure  due  to  any  given  incident  plane  pulse. 
It  is  also  clear  that  the  first  terra  in  Eq.  (9)  represents  the 
reflected  field  due  to  rigid  boundary  condition  and  therefore  the 
second  term  should  correspond  to  the  correction  from  a  rigid  boundary 
condition. 

Assuming  that  the  incident  pressure  is  a  harmonic  pulse  train, 
in  the  time  interval  [0,T], 


Pinc^t)  *  cosuit (H(t)  -H(t-T)] 


(10) 


where  H(t)  is  the  step  function,  the  reflected  pressure  due  to 
normal  incidence  could  be  expressed  as: 


Pref(t)  =  /  Pinc(T>  h(t-t)dT. 
o 


(11) 


Substituting  Eqs.  (9)  and  (10)  into  (11)  one  obtains: 
for  t  <  T 


Pref(t)  =  /  cosuit [<5(t-i)-2iu)0  E  qn  eiU°^C  I^r'n]dT 


or: 


Pref(t)  *  cosuit  -  2i  Z 


9ne 


iui0tnn 


”  <2->  -n2 


U)0'  n 


[e  i'J°tr'n(-inncosuit  +  —  sinuit)+innl 

Uln 

(12) 


and  for  t  >  T 


Pref(t)  =  ~2i  £ 


q  eiuoCrin 

— - - , -  [e  UJ°  nn(-inncosuiT+  —  sinuiT)+innJ 

n  (j-)4V  ,J>0 


(13) 


The  factor  eiuotnn  in  ^q.  ( 13)  represents  an  exponential  decay  with 


A  different  type  of  reflection  that  will  be  considered  in  this 


study  is  due  to  an  incident  FM  pulse.  Let 


PinC(c)  =  cos  [ut(  l+ru»0t)  1  [H(t)-H(t-T)  1 


(14) 


In  terms  of  non-dimensional  quantities,  Pinc^c^  can  exPressefl  as 


wt. 


Pinc(t)  =  cos  [uit(  1+3  — =)  ]  [H(t)-H(t-T)  ] 


(15) 


where  r  =  3  v  and  u  =  ~ 
T 


The  frequency  of  the  FM  pulse  as  expressed  by  Eq.  (15)  is  changing 
from  u)  to  (l+3)w  as  time  changes  from  0  to  T.  Substituting  Eqs.  (15) 
and  (9)  into  (11),  the  reflected  pulse  can  be  expressed  as: 
for  t  <  T 


^ref^c^  “  ^>inc(t^  “  2ito0  Z  qn 

n 


I  qn  ei(i}°trin  /  cos  [ut(  l+rm0T)  ]e  LtJor>T  dT 


or 


pref^c^  =  ^inc 


_  2 

_  (w_nn)  !_i 

(t)  -  2i  I  qn  eiuotnn(i+i)  /llT  {e"1”  *  erf  1 


4  '  „  —  rm  ; - 

2ru 


(r-ict  +  w  -  nn)] 


-i 


(— r— ct 


U) 


0t 

nn>]> 

a=0 


(16) 


and  for  t  >  T 
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m£X2  +  (X2~Xj)k2+  (X2~Xi)c2  =  F 

*'  •  •  • 
m  1 X i  +  Cj  Xj  +  k^  X^  =  (X2~Xi)k2  +  (X2_XOc2 

Substituting  Xj  =  xie^'-ot  and  X?  =  X2ei-'c  one  obtains 

v  klTk2 

X*  =  -  (njiu  +  cj  +  c2  f  — —  ) 


(18-a) 


(18-b) 


where  D  is  the  determinant  of  Fqs.  ( 1 8 -  1 ) 
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D  = 


k2 

(n^iu  +  Q-2  +  ) 

k2 

-<C2  +  -r—  ) 

*■  ia 


k2 

(Iw  +  c2> 

(mjiu)  +  cj  +  C2  + 


The  impedance  function  Z(w)  can  then  be  expressed  as: 


Z(w)  = 


X  j  _  _ 

2  m^io)  +  cL  +  c2  +  ^  (kj+kz) 


Expanding  the  determinant  D  and  dividing  Z(w)  by  pc, 
obtains : 


Z(o>)  =  P(ia) 
pc  Q(w) 

where 


10  i  °°2 

P(u)  -  3|a2(22~  ~  __  _ —  “  w^uj2  +  a^  +  bjb2  +  i 

0)2  ojj  0)^2  ^ 

_  alt>2  a2bl 

aju|b2  -  (— — —  +  - 

wj  a)2 

_  al  a2 

Q(j)  =  (b]+b2)  +  i  [a^yj  -  (-—  +  -_)  J 

al  'J2 


ki+k2 


(19) 


one  finally 


(20) 


[a20)2(bl+b2^  ■*" 


)1 


and 
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Parameters  aj,  a2»  bj,  and  b£  are  dimensionless  positive  numbers.  In 
terms  of  the  new  variables,  the  two  natural  frequencies  of  the 
undamped  oscillating  system  would  be  the  positive  roots  of  the 
following  equation: 

124  -  (w°  +  “2  +  aj  10 1  +  ^“2^  =  °  (20a) 

Obviously  the  shape  of  the  function  | R(u> ) |  would  depend  on  the  values 
of  the  above-mentioned  dimensionless  numbers.  As  it  can  be  seen 
later,  by  proper  selection  of  these  numbers,  we  can  arrange  for 
|R(u)|  to  fall  below  an  optional  level  in  the  frequency  range  of  < 
<*>  <  P-2*  Substituting  the  new  impedance  function  in  Eqs.  (3)  and  (4), 
the  reflected  pressure  due  to  an  incident  plane  pulse  can  be 
expressed  as: 


(22) 
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and 

ala2 

ui  - - 

1  a 

a2 

u2  =  ilal(  l+lj2)+^-  (bi+b2>] 

1  2 

U3  =  (a  +  — )aja2  +  a2  +  ^1^2  +  bj  +  b2 

U4  *  -i[ai(l+b2)  +  aa2(l+b2)] 

U5  =  -aa  ja2 

The  denominator  of  Eq.  (22)  is  a  fourth  degree  polynomial  with 
complex  coefficients  U2  and  U4.  Let  nn  be  the  complex  roots  of  this 
equation  in  the  upper  half  of  the  complex  n-plane. 

u^n1*  +  U2h^  +  U3n^  +  u^n  +  U5  =  0  (23) 

Substituting  Eq.  (22)  into  Eq.  (21),  using  the  Residue  Theorem  and 
closing  the  path  of  integration  along  a  semi-circle  in  the  upper  half 
of  the  complex  plane,  h(t)  can  be  expressed  as: 

<jj°  0 

h(t)  -  5(t)  -  (-—•) ( 2n i )  l  qn  eiUltT1a  (24) 

n 

The  residue  qn  can  be  expressed  as: 

A(nn) 

qn  =  B7!^)  (25) 

and  summation  index  n,  as  mentioned  before,  applies  only  to  those 
poles  located  in  the  upper  half  of  the  complex  plane.  The  complex 
roots  of  Eq.  (23)  can  be  shown  to  be  either  purely  imaginary  or  pairs 
symmetric  with  respect  to  the  imaginary  axis  and  as  a  resuLt  the 
right  hand  side  of  Eq.  (24)  represents  a  real  function.  It  can  be 
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observed  that  with  the  exception  of  the  summation  index  n,  Eq.  (24) 
Looks  exactly  the  same  as  Eq.  (9)  which  was  discussed  before  for  the 
case  of  single  degree  of  freedom  oscillator.  So  we  can  state  that 
Eqs.  (12),  (13),  (16)  and  (17)  would  again  represent  the  reflected 
pressure  due  to  harmonic  pulse  train  or  to  an  FM  pulse  .respectively, 
knowing  that  qn  is  computed  from  Eqs.  (25)  and  (23). 

2.3  Discussion  of  the  Results 


2.3.1  Single  degree  of  freedom  impedance  model 
The  impedance  function  Z(id)  in  this  case  depends  on  two  parameters 
and  bj  expressed  by  Eq.  (5).  The  absolute  value  of  the  reflection 
coefficient  R(u>)  is  plotted  in  Figure  1  for  a  fixed  value  of  bj  and 
three  different  values  of  the  parameter  aj.  The  minimum  value  of  the 
bell-shaped  function  (R(oj)|  depends  on  the  parameter  bi  only, 
occuring  at  the  tuned  frequency  u)  *  u)0  and  its  value  can  be  expressed 
as : 

b  J.- 1 

iR(uo)l  -  ( bprr } 

On  the  other  hand  an  increase  in  the  value  of  aj  would  result  in  a 
sharper  bell-shaped  curve  as  demonstrated  in  Figure  1.  To  show  how 
the  function  R(<j)  influences  the  reflected  pressure,  let's  look  at 
the  following  incident  pulse: 

Pj_nc(t)  =  coswt  (H(t )  -  H(t-T)  J  ,  ^  =  4 

For  a  fixed  R(u),  (aj  *  1.5,  =  1.1),  the  reflected  pressure  is 

plotted  for  different  ratios  of  •  Figure  2  demonstrates  the 

special  case  when  the  frequency  of  the  incident  pulse  train  is  equal 
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to  the  natural  frequency  of  the  impedance  boundary.  Figures  3  and  4 
correspond  to  the  case  when  oj/u)0  =  .8  and  1.2  respectively.  It  is 
observed  that  the  amplitude  of  the  reflected  pressure  could  be 
reduced  to  almost  8%  of  that  of  the  incident  wave  for  u)/uj0  =  1.  Any 
increase  or  decrease  in  the  ratio  of  u)/<jo0  would  lead  to  a  higher 
reflection  amplitude  as  demonstrated  by  Figures  3  and  4. 

On  the  other  hand  for  a  fixed  value  of  u>/uj0,  reducing  the 
parameter  aj,  which  corresponds  to  shifting  to  a  flatter  bell-shaped 
curve,  would  result  in  a  smaller  reflection  amplitude  as  expected. 

In  all  cases,  the  residual  part  of  the  reflected  pulse,  corresponding 
to  time  wt/2ir  >  4.0,  decays  exponentially  as  expressed  in  Eq.  (13). 

Next  let's  look  at  FM  incident  plane  pulse. 

In  Ea.  (15)  let  3  =1,  (J  =0.3  and  u)t/2~  =4.0.  This  represents 
an  FM  pulse  train,  with  a  frequency  range  of  .8  <  u/tu0  <  1.6  for  the 
duration  of  the  pulse.  The  reflected  pressure,  as  demonstrated  in 
Fig.  5,  shows  a  gradual  build-up  in  the  amplitude  with  increasing 
frequency.  Figure  6  represents  the  reflected  pressure  for  the  same 
incident  wave  but  corresponding  to  a  different  value  of  parameter  aj. 
Once  again  it  is  observed  that  a  smoother  reflection  coefficient 
function  would  result  in  a  smaller  reflection.  For  a  fixed 
reflection  coefficient  function,  any  increase  in  the  frequency  range 
of  the  FM  wave  would  result  in  a  higher  reflection  as  demonstrated  by 
Figures  6  and  7. 

2.3.2  Two  degrees  of  Freedom  Impedance  model 

The  reflection  coefficient  as  expressed  by  IP;  s .  (4)  and  (20)  is 
frequency  in  Ficurc  d  for  the  set  of  parameters 


plotted  vs. 


c.cco 


3.CC0 


i 


6.OC0 
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o 

W2 

a  =  —  =  3,  aj  =  1.0,  a2  =  .65,  bj  =  1.0,  and  b2  =  0.7 
*1 

The  natural  frequencies  of  the  undamped  oscillator  obtained  by 
solving  Eq.  (20a)  would  be 

1^2 

—  =  .898  and  —  =  3.338. 

0  0 

“l  “l 

For  the  above  set  of  numbers,  the  absolute  value  of  the  reflection 
coefficient  function  would  fall  below  0.2  in  the  frequency  range 
Qj  <  m  as  shown  in  Figure  8.  The  complex  poles  of  Eq.  (21)  are 

hi  =  (0.0,  .5808) 
n2  =  (.9869,  1.168) 
n3  =  (-.9869,  1.168) 

34  =  (0.0,  6.630) 

An  FM  incident  pulse  represented  by 

pinc(c)  =  cosut(l+  uit)[H(t)  -  H(t-8ir)] 

with  3  =  1  and  iu/u >0  =  1.2  corresponds  to  an  FM  wave  of  constant 
amplitude  and  a  frequency  variation  of  1.2  <  w0 <  2.4 
corresponding  to  the  duration  of  the  pulse.  It  can  be  seen  from 
Figure  8,  the  magnitude  of  the  reflection  coefficient  function  is 
less  than  11%  for  this  frequency  range.  In  Figure  9,  the  reflected 
pressure  has  an  amplitude  of  less  than  20%  for  almost  the  whole 
duration  of  the  pulse.  A  comparison  of  Figures  5  and  9  shows  that  a 
multi-degree  of  freedom  impedance  function  would  minimize  the 
reflection  amplitude  of  an  FM  wave  in  the  frequency  range 
corresponding  to  the  natural  frequencies  of  the  imepedance  model. 


Reflected  Pressure  for  an  incident  FM  Plane 
Pulse  with  .J.  =1.2  and  2=1. 
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Figure  10  shows  the  reflection  for  a  frequency  variation  of 
1  <  (jj/w0  <  3  which  represents  a  gradual  build  up  in  amplitude  to 
almost  45%  around  ut  =  8ir .  Figures  11  through  15  demonstrate  the 
reflection  for  cases  when  the  amplitude  of  the  incident  FM  pulse  is 
not  constant.  Two  types  of  amplitude  variations  have  been 
considered,  linear  and  parabolic.  The  frequency  range  of  the 
incident  wave  is  still  1  <  m/u0  <  3  for  all  these  cases.  The 
reflection  has  been  computed  by  numerical  integration  of  Eq.  (11) 
after  substitution  of  the  residues  from  Eqs.  (24)  and  (25). 
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CHAPTER  III 


PULSE  SCATTERING  FROM  IMPEDANCE  COVERED  CYLINDERS 


3. 1  Statement  of  the  Problem 

In  this  chapter,  the  problem  of  a  pulsed  plane  wave  scattered  by 
an  infinite  cylinder  will  be  solved.  The  cylinder  is  infinite  in 
extent  and  the  incident  plane  wave  is  parallel  to  its  axis.  The 
cylinder  is  assumed  to  have  a  cover  that  possesses  a  locally  reacting 
impedance  as  the  case  in  Chapter  II. 

Consider  a  plane  pulse  traveling  with  speed  c  in  the  x-direction 
perpendicular  to  the  cylinder's  axis  (see  Figure  16). 

Pine  -  Po<5(f  -  t)  (26) 

Applying  complex  Fourier  transform  on  P^nc,  one  obtains: 

Piu  -  I  Pine  e1UC  dt  *  /  P0  6(|  -  t)eiut  dt  =  PQ  elkx  (27) 


— oo  —oo 

where  k  =  rn/c  is  the  wave  number.  The  incident  wave  can  be 
reconstructed  by  applying  inverse  Fourier  transform: 


t) 

(28) 


The  integrand  of  Eq.  (28)  Pj_u  e~iwt  represents  a  time  harmonic 
component  of  the  incident  wave.  Since  the  scattered  field  due  to 
this  component  of  the  incident  wave  would  have  the  same  time 
dependency,  factor  e“lu)t  for  convenience  will  be  suppressed.  P^w 
can  be  expanded  in  terns  of  cylindrical  functions  [17]  as: 


=  Pr 


ikreos? 


where 


PQ  E  en  in  Jn(kr)cosno 
n=0 


(29) 


z 


0  = 


1  and  t 


2  for  n  *  0 
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The  scattered  field  PS0J  corresponding  to  an  incident  pressure  P^w  can 
be  expressed  as 

CO 

P su>  =  Po  £  An  en  in  H^(kr)cosn$  (30) 

n=0  n 

the  Fourier  coefficients  ^  are  obtained  by  applying  the  boundary 
condition 


piw  +  psw 
uiw  +  usm 


=  Z(w)  at  r  =  a 


where  u^w  and  uS(l)  represent  the  incident  and  scattered  components  of 
velocity  in  the  radial  direction  and  Z(ui)  is  the  impedance  of  the 
boundary,  which  is  a  function  of  the  frequency. 

Ui^  and  usw  can  be  expressed  as  follows: 


1  3  piw 


P  00 

ro 


imp  3r  2pc  _ 
n=0 


I  en  i  [J  (kr)-J  (kr) jcosnp 
„  n  n+1  n-1 


1  3psm 

Usa  =  ~ 


p  00 

~2~~  £  An  en  in+1  [H^|(kr )  -  H^Ckr )  ]  cosnp 

n=0  (32) 


where  p  is  the  mass  density  of  the  acoustic  medium  surrounding  the 
cylinder.  Substituting  P^,  Psu>  u^w  and  usu  in  fc)q.  (31)  and  solving 
for  the  constant  An,  one  obtains: 


An(ka) 


Jn(ka)  +  j'(ka) 

11  pc  n 

H(  °(ka)  +  -Z(  -  H'U)(ka) 
n  pen 


n  =  0,  1,  2  ...  (33) 
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where  the  prime  stands  for  differentiation  with  respect  to  the 
argument  of  the  cylindrical  functions.  Finally  the  scattered 
pressure  due  to  an  incident  plane  pulse  on  the  cylinder  can  be 
expressed  as: 

.  +°° 

PsCr,t,4>)  PSUi  e“1Ut  dm, 


P  00 

rO 


*  jr  E  Sn^08^  /  An(ka)H^ 1  ^ (kr)e~iu)tdu). 


Letting  n  -  ka  and  £  -  c/a,  then  Eq.  (34)  can  be  written  as: 


Pq£  00  -fCO 

Ps(r,t,<>)  =  E  enincosn  <pf  An(n)H(  U  (n  -  )e“irKCdn  (35) 


n  a 


3.2  Impedance  Function 

Consider  the  simple  oscillator  as  a  model  for  the  impedance  of 
the  surface  expressed  as: 


m  0  2  k 

Z  Cq  imw0(—  )  ,  cm  =  — ■ 

°  0  m0  m  o  m 


Dividing  Z  by  pc  and  introducing  the  constants 


bi  =*  —  and  a>  =  /—  , 

1  pc  1  km/ pc 


*-*  t  .  /  LJ  \ 

—  =  bj  -  iai( — - ) 

PC  *  1  m0  m 


To  express  Che  impedance  in  non-dimensional  variables ,  the  following 


variables  are  introduced: 


_  2ttc  J  =_  = 

“2.1 t  a 


Then  the  variable  n  can  be  written  as 


n  =  ka  =  (—)(—) 


w°  I 


or  (— )  =  nX 


Finally  let  Z(n)  = 


Z(n)  =  bj  -  iaj(An  -  — 


The  impedance  function,  as  expressed  by  Eq.  (36),  is  expressed  as  a 
function  of  the  radius  of  the  cylinder  as  well  as  the  characteristics 
of  the  oscillator.  Finally  the  reflection  amplitude  coefficients  are 


given  by: 


j  (n)  +  iZ(n)J  (n) 
11  n 


An(n)  =  - 


n  =  0,  1,  2,  ...  (37) 


H(1)(n)  +  iZ(n)H  (1)(n) 
n  n 


Using  the  identify  Jn(n)  =  1/2  iH^^n)  +  Hn ^ ^ ( ri )  ]  ,  An(':)  can  be 


rewritten  as: 


H(2)(n)  +  iZ(n)H,(2)(n) 
An( n)  -  (1  +-2 - - - 

H(  L)(n)  +  iZ(r.)H,(  i}(n) 
n  n 


-  7  11  +  Bn(n) 
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<v 

. 


k./ 


t* 


P 


►« 


Substituting  Eq.  (38)  into  Eq .  (33)  one  obtains: 


Ps(r  ,t  ,$)  =  —7 -  I  en  incosn?  /  —  )e  ^^dn 

H'n  ,,  “  J  n  a 


n=0 


P  Z  00  _^CO 

-  4^—  ^  enincosn.?  /  Bn(n)H^(n  •-  )e  instdn  (39) 

n=0 


It  can  be  shown  that  the  first  integral  on  the  right  hand  side  of  Eq. 
(39)  has  no  contribution  to  the  scattered  pressure. 

Let 


Q  =  Z  en  incosnj>  /  H^^(n^-) 

14  '  n  a 

n=0 


.  .  ct 

+“  , .  x  -in  — 

a  j 

e  dn 

n  a 


first  let's  consider  the  case  when  r/a  »  1.  Using  the  asymptotic 
expression  for  H^^nr/a)  one  obtains: 

„  r  .n  . ,  /  2~  x  -i(n7i/2  +  n/4)  ^  1  -in-—-- 

Q  ~  Z  eni  cosno(/ - r-  )  e  J  -  e  a  dn 

x  „  n  f  irr/a  ‘ 


n=0 


/n 


or 


A  y  2  -i7r/4(-+00  1 

Q  =  /  -  e  I  —  e 

711:73  -  K 


but 


.  ct-r 
-in- 


—  ,  /  x  .n  -iniT/2, 

a  dii  i  Un  cosn$  1  e  } 

n=0 


Z  £n  cosnoi  e 


.n  -in:r/2  „  „r  1  ,  .  , 

‘  ~  =  t.  tn  cosno  *  2{  -r  +  cosd>  +  cos2$.  +  ...} 

n=0  n=0 

=  24  -  4^  =  0 

L 


hence  Q  =  0.  Furthermore,  one  can  consider  the  integral 


I „-/  .‘‘"dn 


(u  =  ct/r  >  1) 
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This  integral  is  obtained  by  a  change  of  variable  in  the  expression 


for  Q 


Q  =  -7—  Z  cn  in  cosn^(In) 
r/a  n=0 


it  can  be  proven  that  Q  is  a  pure  imaginary  number  and  hence  would 
not  have  any  real  contributions  to  the  scattered  pressure.  Let's 
consider  the  following  contour  integration 


Jc  <1}(n)  e-ictndn , 


where  C  is  the  closed  contour  shown  in  Figure  17.  The  function 
Hn^(n)  does  not  have  any  pole  inside  the  closed  contour. 
Furthermore,  on  V  the  integrand  is 

H(I)(n)  eIc‘n  -  Jlht  e+i<!’  -  T  *  e~1“n 

n 

-  STTZZ  ~  r)  -m(o-l) 

=  /2/7rn  e  2  4  e 

since  (a-1)  >  0  then  e-in^a-^  vanishes  on  F  .  On  the  small  contour 
cz >  |h|  «  1«  The  asymptotic  expressions  for  the  Hankel  function  for 
small  argument  [18]  are  given  by: 

u ( 1 ) /  \  2i  „ 

H  (n)  ~  —  *nn 
o  u 

3nd  H(1)(n)  ~  -  -  -(-n—  ^  2*  for  n  >  1 

^  ^  _  n 
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Consider  the  n  =  0  terra*  letting  n  =  ee^ ,  then 


Io  -  /  H^1)(n)e"ia%  =  -  '  H<;u(n) 


(1),  N  -ian  . 
(n)e  an 

c  o 
e 


•  i9  ,  ° 

-  -  J  —  in(eeid)  e"laCe  (ieeib)d9  =  —  J  (10  +  J£ns)ei0de 


~7T 


4i 

- - eine  -  2e  =  0  as  e  -*•  0  • 

77 


for  n  >  1  the  integral  is  given  by: 

In  =  I  H^1)(n)e'iandn  =  -  J  (H<0(n)e 
-00  £ 


-lan 


]dn 


|  ( 2n)  ((-^{|!-  e 1-n  [1  +  (-l)nl  - 


0 


n  =  odd 


_(i)2n+l(n-2)!  el-n 
n  «  even 


Summing  up  all  these  contributions,  then 


Q  =  —  Z  en  in  cosn?  In  =  (—)(—)  E  (real  fun.) 
r  n=0  r  71  n 

is  purely  an  imaginary  number. 


3.3  Watson  Transformation 

The  scattered  pressure  due  to  an  incident  plane  pulse  on  the 
cylinder  can  thus  be  expressed  as  : 


Pc(r  ,t  = 


+» 


n=0 


cosnj  j  Bn(n)H^^(n  — )e  inct//a  dn 
_  n  a 
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where 


H^2)(n)  +  iZ(n)  H^(2)(n) 

n  H(1)(n)  +  11(11)  H,(1)(n) 

n  n 


Let  En(n)  =  Vn)  H*1}(n 


Using  the  following  properties  of  the  cylindrical  functions 


Hl»(n)  -  ^  n‘U<r,>  , 
H*2,(n)  -  e"Y’i  H<z)(n), 


it  is  easy  to  show  that  E_y(h)  =  e-Y^i  Ey(n).  For  the  special  case 
that  7  is  an  integer  one  obtains: 


E  (n)  =  (-Dn  E  (h) 
-n  n 


Eq.  (41)  would  lead  to  the  following  identity 


I  en  i  cosn$  En(n)  *  £  incosn$  En(n) 

n=0  n= 


Using  Sommerf eld-Watson  transformations,  the  right  hand  side  of  Eq. 
(42)  can  be  expressed  as: 


l  incosn<?  En(n)  =  \ 


j-yfi/ 2  C0Sy;  ^(n) 


sinyTT 
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where 


,d  sinyir. 

^  dy  ; 


(-lA 


Y  =  IT 


with  the  contour  C  shown  in  Figure  18.  It  can  be  shown  that  the 
value  of  the  integral  on  the  lower  half  of  the  contour  'C'  is  the 
same  as  that  for  the  upper  half.  Letting 


I  -  / 


-yiri  ,  _  ,  , 
e  cosy<!>  Ey(n) 

C  sinyir 


dY  =  /  (  )dy  +  /  (  )dy  (44) 


and  making  a  change  of  variable  (y  +  -y),  one  obtains  an  expression 
for  I  as: 


+<*> 

I  -/  ( 


)dy  +  /  cosy<>  E_y(o) 


sinyir 


which,  upon  substitution  of  E_Y(n)  =  e-^1^  Ey(n)  becomes; 
+“  -yirj 

I  ~  2  f  e  2  cosytt  E  (n)  — ~ — 

J  Y  smyTT 


Incorporating  the  above  transformations  into  Eq.  (40),  the  scattered 
pressure  field  could  be  written  as: 


„  /  ,  ,  .  P°^  7”  cosy;  Sv(n)  ,,(1),  r,  _<IU,r/3J  : 

Ps(  r  ,  t ,  y )  -  -1  -r-  J  /  - - - - -  Hy  (n-)e  mct/adYdn 


Y  a 


where 


( i)  —  '  (?) 

H^;(n)  +  iZ(n)  A) 
BY(n)  =  “  •  (" n 

H.j, 1 ; ( n )  +  iZ(n)  H ^ '  w(n) 


(47) 
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Following  Uberall  [7J,  the  scattered  pressure  as  expressed  by 
Eq.  (46)  could  be  split  into  two  parts: 

1)  A  geometrically  reflected  field 

2)  Creeping  or  circumferential  waves  that  encircle  the  cylinder 
any  number  of  times  and  propagate  to  the  observer. 

The  presence  of  circumferential  waves  has  been  confirmed 
experimentally  by  measuring  the  multiple  echoes  that  return  in  the 
scattering  of  a  single  pulse.  This  can  be  demonstrated 
mathematically  by  expanding 

cosy<p  =  e  cos{(it  -  $)y]  -i  e  sinyTr 


1 

sinyn 


iYTr  r 
-  2i  e  l 


2isy^ 


The  product 


e-yni/2  cosYd/sinYTT 

can  thus  be  expressed  in  either  of  the  following  forms: 


and 


or  in.  the  form: 


-Yiri 


m 


e  cosyft 


sinyTr 


-  [-i  e  *  Z 

\  =  ± 


Z  e 
s=0 


+  2tts)j  _ieiy(-2  -  <t>) 


(48-b) 


where  <j>  =  2tt  —  <J>  , 


and 


<J>+  =  4> 


Substituting  Eq.  (48)  into  (46)  one  obtains: 


Po? 


ps(r  »t  ,<))  =  -  —  (Pj  +  PG) 


(49) 


where 


pUc.c.t)  -I  C  r  r  e1T(^  +  2’s\(n)n‘u(n 

X=  ±  s=0  —  -»  ’ 


:/ a 

c 

(50) 


and 
G 


p  (r ,t ,$)  «  J  /  e1Y(2  ~  |)aY(n)H^1)(n  |)e  inct/a  dydn  .  (51) 

— oo  —co 


In  Eq.  (49),  corresponds  to  geometric  optic  reflection  and  p^c 

corresponds  to  creeping  waves  on  the  surface  or  close  to  the  surface 

of  the  cylinder.  The  angle  $  ’  represents  the  waves  that  encircle 

+ 

the  cylinder  in  c.c.w  direction  and  the  angle  $ '  corresponds  to 
c.w  creeping  wave  (see  Fig.  19).  These  waves  could  be  observed  both 
in  the  insonified  and  shadow  boundary  of  the  cylinder. 
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The  summation  on  s  indicates  the  presence  of  creeping  waves  that 
encircle  the  cylinder  more  than  2tt  radians  in  the  positive  or  the 
negative  sense. 

On  the  other  hand,  substituting  Eq.  (48)  into  (46),  the 
scattered  pressure  can  be  written  in  a  different  form  which  would  be 
valid  for  an  observer  in  the  far  field: 


P°5  c  G 

Ps(r  ,t  ,4>)  =  -  <Pn+  P  ) 


(52) 


where 


'II 


A=  ±  s=0 


os  +oo  +CO  >  I 

I  S  1  *  2  *lY(^  +  2"\(n)Ha,(r,  — )e-inct/adtdn 


(53) 


and  p^  is  again  given  by  Eq.  (51).  In  Eq.  (53),  pc  represents 

II 

circumferential  waves  travelling  to  an  observer  in  the  farfield  in 
either  the  insonified  or  the  shadow  region  (see  Fig.  20). 

If  t  =  0  corresponds  to  the  arrival  of  the  wave  at  the  shadow 
boundary  then  the  arrival  time  for  an  abserver  at  r/a  »  1  and  angle 
$  (see  Fig.  21)  could  be  expressed  as: 


r  =  J a?-  +  r2  - 


2arcosu 


IT  - 


r  3  r/l  +  (— ) 


2a  .  0 

—  sin 
r  2 


=  r  -  asm— 


Figure  20a.  Circumferential  Waves  Traveling  to  a  Far 
Distant  Observer  in  1$)  >  tt/2  Region. 


Figure  20b.  Circumferential  Waves  Traveling  to  a  Far 
Distant  Observer  in  <  ~/2  Region. 


Figure  21.  Plot  of  the  Geometrically  Reflected  Pressure  for 
an  Observer  in  the  Far  Field. 


Then  Che  arrival  time  for  the  reflected  pulse  can  be  expressed  as 


5(ct  -  r  +  a  sin  y)  =  6(ct  -  r  +  2a  sin  -|o 

For  a  >  45°,  the  geometrically  reflected  wave  reaches  into  the  (J  < 
90°  region. 

3.4  Geometrically  Reflected  Pressure  (r/a  »  1,  0  <  <{>  <  n) 

Since  the  reflected  pressure  is  symmetric  with  respect  to  the 
x-axis,  angle  <p  could  be  limited  to  0  <  £  <  .  Dropping  the  factor 

-P0^/4tt  for  the  time  being,  the  geometrically  reflected  pressure  is 
given  by 


pC(r,t,})  -  r  r  «lY<^  "  *)BT(n)«i1)(n  f)e-in£S  dYdn 


(54) 


where  B^(n)  is  expressed  by  Eq.  (47).  Letting  n  »  y  secct,  then  the 
cylindrical  functions  Hy^^(n)  and  Hy^^(h)  can  be  expressed  in  the 
following  generalized  form: 

H^°(Yseca)  =  FX(y .cOe1'^  ’a) 

and 

H.’(l)(Yseca)  =  F2(Y,cO  'a)  (55) 


*(Y,a) 


Y(tana-o)- 


TT 

T  • 


where 
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For  the  case  when  |y|  »  1,  Fj(y,a)  and  F2(y,a)  can  be  expressed  as 
(18): 

Fi(y,a)  =  /  2/ Try  tana  (L(y,a)  -  iM(y,a)) 

F2(y.a)  =  / sin2a/7rY  [iN(y,a)  -  0(y,a)] 

where  L,  M,  N  and  0  are  expressed  in  terms  of  the  following 
convergent  series: 


L(y ,a) 


Z 

k=0 


U2k( icota) 

2k 

Y 


M(y ,a) 


-i  l 
k=Q 


u2k+l(icota) 

^2k+l 


N(y,a)  s  z 
k=0 


^2k( icota) 
Y2k 


0(y ,a) 


i  I 
k=0 


V 

2k+I(icota) 

2k+i 

Y 
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and  uQ(t)  =  1 

Uk+l(t)  =  I  t2(1  "  t2)uk(C)  +  8  n-5t2)uk(t)dt  k=0il)2>>>> 
VQ(t)  =  1 

vk(t)  =  uk(t)  +  t2(t2-l){  j-  uk_^(t)  +  tuk_1(t)}  ,  k=l ,2 , . . . 

Substituting  for  y  in  terms  of  the  variable  a  and  n  (Y=ncosu)  one 
obtains : 


H^2)(n)  +  iZ(n)H,'(2)(n) 
ncosa(n)  H(D(n)  +  iz(n)H^(1)(n) 


Y  =  ncosct 


F  (ncosa ,a)  +  iZ(p )F  (ncosa ,a) 

1  2  -2i0(ncosa ,a) 

- -  e 

F^ncoso^a)  +  iZ(n)F2(ncosa ,a) 


=  G(a,r\)e 


— 2  i 


(55) 


where  +>(ncosa,a)  =  n(sina  -  acosa)  -  tt/4. 
For  an  observer  in  the  farfield  (r/a  »  i) 


u(l),  r,  j  2  i  ( n  r  /  a - j  —) 

H  (n  — )  ~  / — —  e  2  >4 

Y  a  it  r/a 


Eq.  (54)  can  now  be  expressed  in  terms  of  the  new  variable  as: 


56 


'(r  ,t  ,9)  -  /  /  C(a,n)e~2iil'  e,lncos“<T  '  V  .  f~2 


t\=  — 00  '-a 


unr/a 


.  /  L  TT  It  .  .  Ct 

.  lCn  —  -  ■=■  n  cosa  -  7*)  -in — 

•  e  a  2  &  •  P  a  •  f-nci 


4  •  e  a  (-nsinadot)  dn 


pG  (r,t,<>)  -  -  /  —  ^  e**^  /  /n  e  ^n"dn  /  G(a,n)sinae  ir'^^a^da 

tt  —  n=  C  a 


where 


g(«)  =  2sina  +  (<+>—  2a)cosa  , 


and  Ca  describes  Che  path  of  integration  in  complex  ct-plane  (see  Fig. 
22).  Changing  the  direction  of  integration  Eq.  (56)  can  be  rewritten 


iTT  +■» 


P  (r,t,$)  =  e  ^  ^  Kh)  e  in“dn 


where  I(n)  =  /  »  G(a,n)sina  e  in®^a^da 
Ca 


Eq.  (59)  can  be  approximated  by  integration  along  t.  e  path  of  the 
steepest  descent. 


3.5  Method  of  the  Steepest  Lescent 

The  saddle  point  for  the  integral  in  Eq.  (59)  is  located  at 
a0  where 


g’(a0)  =  0 


or  ~  =  (2a  -  <p)  sina  =  0,  a0  =  4 
a  a  J  l 


The  equation  of  the  steepest  descent  path  (Fig.  23)  is  expressed  by 


Im[g(a) ] 


const  =  Im[g(a0) ] 


Im[2siny]  =  0 


Substituting  a  =  ar  +  ia^  ,then  the  equation  becomes: 

( 1  -  cotha^)  *  (y  -  ar)tanar  (60) 

The  asymptotes  of  the  steepest  descent  path  given  by  Eq.  (60)  are 
a^  -*■  ±  00  at  ar  =  tt/2. 

The  slope  of  the  steepest  descent  path  at  a  =  a0  is  given  by: 


da^ 

dar 


da^ 

^dar 


-tanar  +  (-|-  -  ar)sec2ar 


-cotha^  +  o^/sinh2^ 
-tana/2 


a  =  a0 


It  can  also  be  shown  that  the  integrand  in  Eq.  (59)  vanishes  on  Fj 
and  T2*  Thus,  the  integral  in  Eq.  (59)  is  seen  to  be  equivalent  to 

I(n)  =  /SDpG(a,ri)sina  e  inS^da  (61) 

Letting  W2  =  g(a0)  -  g(a)  , 

or  W2  =  2sina0  -  2sina  +  2(a  -  a0)cosa  , 
dW 

then  2W  —  =  -2  cosa  +  2cosa  -  2(a  -  n0)sina  =  -2 (a  -  a0)sina  , 

,  da  W 

and  ~rr  ~  ~  7 - : — : - 

dW  (a-a0)sina 


Substituting  for  a  in  terms  of  W,  Eq .  (61)  could  be  expressed  as: 


Saddle 

Point 


Ficure  23.  Steepest  Descent 
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w 

1(a)  -  -2  /  iG<a.n)W(q) j  g-in [g(a0)_W? ] 
W=0  a  ~  aO  a=a(y)  C 


=  2e~^^Tls^'ncio  r  ^G(q,n)W(a) .  ginW2 

r  t _ r\  . 


a=a(W) 


Letting  F(W,n)  =  r G( ^ ^ )Vv'( a) . 

Ct  —  J 


a=a(W) 


then  F(W,n)  could  be  expanded  in  a  Taylor's  series  about  W  =  0(a=ao) 


as  follows: 


F(w,n)  -  MO.,)  +  |ij  F(w,,)j  »  +  (£  F(IM„  «i  +  ... 

f  r_rv  .0  /  • 


with  ~  F(W,n)  =  -i-  f-G-(gt.»n)W(a)  ■  dot 
dw  dct  1  a  -  ao  ‘  dW 

-  {G'(a,n>^  +  G(a  }  d_  rj£(2))}  da 

a-ao  ,n'  da  la-a0ji  dW  * 

On  the  other  hand  expanding  g(a)  in  a  Taylor’s  series 


W2  =  g(a0)  -  g(a) 


(a-aD)2 


~  S(«0)  ~  (g(a0)  +  (a-a0)g '  (a0)  +  — -  -  -  g"(a0)  +  ... 


“  <a~ao)2  l~  2T~g’ ’ (a0)  -  4to ' '  '  (a 


0  ^  ^  ^*”^0  )  •  *  •  J 


^a-a0^  “  ^~2!  ®  ’(ao^  ~  37  '  '  (ao^  (a-a0)  • . . ) 


which  is  approximated  in  the  limit  when  W  -  0  (a  *  a0)  by: 


i  •  /  w  ,  ,2  1/9  2sina~ 

llra(  )  .  (-^g.'Uo))1'2  „  =  /— 


sma0 
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Differentiating  Eq.  (63)  with  respect  to  a  one  obtains: 


d  1  I#/  >  2  /  xXV/  x  if  1  it/  \ 

-  2  I_  IT  s  (cto)  "  TT  («“ao)g  (ao>  +  11-  —  s  <ao> 


^y(a-a0)g"'(a0)  +  ...J  1/2 


, .  d  ,W(a) s 

or  lim  —  ( - ) 

da  a-a0 

a  *  an 


j  I"  l  g,,,(a0)J[-  j  g”(a0)] 


3/-sina0 


Differentiating  Eq.  (63)  once  more,  then 


lira  —  (-^— )  =  - 

da2  «-«o  _ 


°  ,1^1  ?  \ 
( 4  +  9  cot^a0) 


a  -*■  ar 


Finally  an  expression  for  first  derivatives  of  F  is  given  by: 

d  , _  -cosa0  _j  _ 

[—  f(W,n)l  =  [G'(a0,n)  /-sina0  +  Gta^.n)  -  (fir — )  •'-sinao 

dW  W-0  3/-sina0  Sin“° 


=  G'(a0,n)  + -j  cota0  G(a0,n) 


Similarly 


F(W,n)  = 


d  rr/»,  „ \ W ( a )  ,  _Nd  /W(a)N1/daN1  da 

_  .  [G  (a,n)—  +  G(a,n)^  (— JJ^))  517 


=  { [G' '  (a,n)  +  G'  (a ,n  )  ~~  (~— )  +  G'(a,n)~-(-^-~) 

a-a0  da  a-ao  da  a-a0 


s  d2  ,W(a) . ,  da 
+  G(a,n)  -  (~-  ---  )  1  -jrr 

da2  a_a°  d,< 


^F(W.n)  ,  , 

/Q» _ v  d__  /day  da 

^  da  ;  da  dW 
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and 

d  ,dcu  d__  ,  W(a)  ,  _  _  d  , W(a)  ,  1  +  W(a)  cosa 

da  ' dbr  da  (a-a0)sina  da  a-a0  sina  a-a0  sin2a 


or  [- —  F(W,n)]  =  {G' ' (a0 ,n)  /-sina0 
dW2  W=0 


2G '  (a0  , T) ) 


-cosa 

3/-sina0 


sina0 

+  G(a0,n)[ -  ■  - 

/-sina0 


<4  +  4  cot2ao>)K- 


/-sina0 

sina0 


2 

) 


cosar 


+  [G'(a0,n)  +  -j  cota0  G(a0,n) j l 


(- 


- -  sma0 

j/-sina„  u 


)  +  /-sina0 


cosa0 


sin2a0 


Finally 

l—  F(W,n)  J  =  {G”(a0,n)  -  (t  +4  ^oc2cto)G(a0>n)}C  -~1- - )  (66) 

dW2  W=0  /-sina0 


The  Taylor's  expansion  of  F(W,n)  about  W  =  0  could  be  written 
as : 


F(W,n) 


cota0 

G(a0,n)  /-sina0  +  (G'(a0,n)  + — ~ —  G(a0,n)]W 


+  [G"(a0,n)  -  <4  +4  cot2a0)G(a0,n)](  — - - )  |j-  +  ... 

/ -sina0 

(67) 


using  the  relationship 


dy  =  / 2n / p 


(2n)’ 
pnr(n+l ) 


the  integral  over  the  steepest  descent  path  becomes: 
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The  function  G(a0,n)  is  of  the  order  0(n).  In  order  to  apply  the 
reside  theorem  to  the  integral  in  Eq.  (70),  it  is  necessary  to  write 
G(a0,n)  as: 

G(a0,n)  =  -1  +  H(a0,n)  (71) 

where  the  function  H(a0,n)  is  of  order  U(l/n)  for  |n |  >>  1.  From  Eq, 
(55) 

(Fj  +  Fj)  +  iZ(n)(F2  +  F2) 

Fj  +  iZ(n)  F2 


H(a0 ,n) 


(72) 


and  for  the  special  case  |y|»l,  it  can  be  seen  that 


H(a0,n)  =  2 


L(a0,n)  -  iZ(n)0(ao ,n)sina0 

[L(a0,n)  -  iM(aQ ,n )  ]  +  iZ(n) [iN(a0 ,n)  -  0(ao ,n ) ]sina0 

(73) 


Substituting  the  expressions  for  L,M,N  and  0  in  the  above  equation 
and  letting  y  =  ncosa0,  it  can  easily  be  seen  that  H(a0,n)  is  of 
1/n.  Pulling  out  the  factor  (-1)  from  G(a0,n)  corresponds  to 
separating  the  major  contribution  due  to  rigid  boundary  condition 
because : 


“2i<> 


•  G(a0,n)  =  [-7 

t  A 


H^2)(n)  +  iZ(n)H^(2)(n) 


H^A)(n)  +  iZ(n)H^(i)(n) 


<<2)(n) 

»;(1)(n) 


1  + 


<(2,(n) 

n;(1)(n) 


-w(H^1)(n),il^2)(n)l  H^(2)(n) 

(H^I}(n)  +  iZ(n)H’(1)(n)]H^(1)(n)  H^(i;(n) 

_ 4i/Tin _ 

lH^°(n)  +  iz(n)H^(1)(n)]H^(1)(n) 

where  W  is  the  Wronskian  and  y  =  ncosa0.  The  first  term  on  the 
right  is  of  order  1/n 2  for  j  n  j  >> 1 .  The  second  terra  corresponds  to 
reflection  from  rigid  cylinder  and  can  be  rewritten  as: 


<(2>(n) 

<U,(n> 
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0(ao,n)  +  iN(a0,n) 

0(ao,n)  -  iN(a0,n) 

-2i'ip(a0  ,n) 
e 


=  [-1  +  2 


0(ao ,n) 

U(a0  ,n)  ~  iN(a0  ,n) 


-2iy(a0  ,n) 
e 


The  second  term  in  the  bracket  can  be  recognized  to  be  of  the  order 
( 1/n)  and  hence  the  term  (-1)  represents  the  major  contribution  to 
the  reflected  pressure  due  to  a  rigid  cylinder.  Finally  from  Eqs. 
(70)  and  (71)  one  can  write  an  expression  for  the  geometric  optic 
terms  as: 


pG(r,t,$) 


[-1  +  H(a0,n)]e"ir,(^  +  2sina°) 


ct — r 

+00  -in( -  +  2sina0)  +*> 

{  /  e  3  dn  -  /  H(a0,n) 

-.00  —CO 


-in(H  +  2sina0) 

e  dn  (74) 


The  first  integral  can  be  recognized  as  a  delta  Dirac  function.  The 
second  integral  can  be  computed  for  (!2  +  2sina0)  >  0  by  closing  the 
path  of  integration  along  a  serai-circle  of  radius  infinity  in  the 
lower  half  of  the  complex  n-plane.  The  poles  of  the  integrand  are 
the  roots  of  the  denominator  in  Eq.  (72)  which  are  the  same  as  the 
roots  of  transcendental  equation: 
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H^(n)  +  i2( n ) ^  (n )  =  o,  for  Y  =  ncosa0  (75) 

Let  nn  be  the  complex  roots  of  Eq.  (75)  in  the  lower  half  plane. 

Then  Eq.  (74)  can  be  written  as: 

~inn 

.  ~  .  r.  «  - (2sin<}>/2+ct-r) 

ct-r)  +  2m  2  2.n  e  a  Y 

n 

.  u(£il£  +  2sin4>/2)]  (76) 

a 

where  u(t)  is  a  unit  step  function  and  £n  are  defined  by 
(F^  +  F|)  +  iZ(n)(F2  +  F2) 

^n  =  —  at  n  =  hn  (77) 

^  (Fl  +  iZ(n)F2| 

To  express  the  residue  £n  in  terms  of  cylindrical  functions  one  may 
write 


P  (r,t,$) 


/2sin6/2U5(2aSini 

r/a 


Since  Che  order  of  Che  Hankel  funccion  is  dependenC  on  n,  Chen  Che 


dif f erenCiaCion  should  be  applied  to  both  the  argument  and  Che  order. 
Let  a  prime  represents  dif f rentiation  with  respect  to  the  argument  and 
a  dot  stands  for  differentiation  with  respect  to  the  order.  Using 
the  chain  rule  of  differentiation,  then: 


R  =  {h'(1)  +  iZ  H  +  iZ  H  +  (h'(1)  +  iZ  }e 

Y  Y  Y  Y  Y  dn 


n  =  n 


n 


=  {(1  +  iZ)H^(1)  +  iZ  h”(1)  +  4-  iZ  H' 1))cosa0}  e  ll^(ao>rin) 

{  T  Y  n  =  nn 


Using  the  identities: 


H''(1)(n)  =  -  rH*(1in)  +  (—  -DK^Cn)  =  - 


n  Y 


2  Y 


^  -  sin2a0H^1)(n) 


or : 


H’(1)(n)  =  (n)  -  ^  H(1)(n) 


H*’(1)(n)  =  -^H^U(n)  H^(1)(n) 


then  R  becomes : 


— )h'(1)-  -■•7- 
n  y 


R  =  {(1  +  iZ  -  —  )H  vw-  iZsin2a0H  ^  cosa0  [(1  - 


I  •  iZ(n))H^(1) 


+  iZ( n ) ( H*  J ^  ^  H^)]}  e  lv(ao»hn) 


Substituting  for  H  (nn)  =  “  - - - H  (nn)  for  7  = 

Y  iZ(nn)  Y 

then  R  can  be  expressed  by: 


R 


~7  1  —  uut>u0  ,  .  . 

{ [— — - +  —  -  iZ(sin2ct0  +  -  )  ]  ir  '  (nn)  +  cosa0 

Z(nn)  Z  "n  ""  Y 


n'<l)(n„)  +  iZ(nn)cosa0  h'*}*]  •n"> 


or 

R  =  {T(a0  ,nn)H^(nn)  +  [l-iZ(nn)c°scto  Jcosa0  H  ^(nn) 


+  iZ(nn)cosa0  H.(^|\nn)}  *  e 

Y  =  ncosa0 

where 

1  _  cosao  j 

T(a0,nn)  - -  [i  -  z  (n0)  -  iZ(nn) [sin2a0  +  — — -J  +- 

z(n„)  nn  n’ 


From  Eq.  (77),  the  residue  &n  can  be  expressed  as: 


ln  =  w - — UJ  (n)cosi(/  +  Y  (n)sinij/]  +  iZ(n)[J  (n)cos^  + 

*  r*n '  ’  I  i 


nncosa0 


( l-icosa0Z) 


(78) 


Y  (n)sin’;)} 
n"nn 
(79) 


The  derivati%re  of  the  Hankel  function  with  respect  to  its  order 
which  is  needed  in  computation  of  Eq.  (78)  can  be  obtained  from: 
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~  H^(x)  =  H  ^(x)  =  (( l+cotTy)£a(^) - - -  ]J  (x) 

T  T  sin2Try  ' 


+  __i -  [£n(^)  +  TTcotY")  J  v  (x) 

sinYx  2  ~Y 

2  k 

wx,  *  ,.k  «(Y+k+l)  l4  ' 

(  1+COtYTT  )  (2  )  2  (  1)  r(  +  k  +  k! 

k=0 


x2  ^ 

-  j0  <-»k  «<» 


The  functions  i[i  in  Eq.  (80)  could  be  obtained  from  its  asymptotic 
representation  [18]: 

00  Bo„ 

(81) 


(,(Z)  .  (*nz  -y  -  : 


n=l  2nZ‘ 


where  Bn  are  Bernouli's  numbers  [18].  The  following  asymptotic 
series  for  y(Z)  can  be  used  for  large  values  of  |Z|: 


V(Z)  =  UnZ  -  (  — — 

22  12Z2 


120Z4 


1 

252Z6 


1 

240Z3 


+ 


...) 


The  following  relations  could  then  be  used  to  compute  Z)  for  any 
values  of  Z 


*U-D  =  v-(z)  -  ytt 

V(-Z)  =••  ~(Z)  +  Y  +  ^TCOtTZ 
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Finally,  multiplying  by  the  factor  -P05/4ir,  which  was  dropped  out 
earlier,  one  obtains: 


pG(r ,t,p) 


PD  / 2sind/ 2  , 

/  - - -  {cS(2asin4  +  ct-r) 

4n  r/a  2 


+  2iri(— )  Z  Zn 
a  11 
n 


(2asin^-  +  ct-r) 


u(f2  +  2sin-^-)} 


(82) 


f:ven  an  incident  plane  wave  Fin(t),  one  may  apply  Duhamel's  integral 
on  Eq.  (82)  to  compute  the  geometrically  reflected  fic.  d  for  any  given 
incident  plane  pulse.  For  a  cosine  pulse 


Finc(t)  =  cosmt [ H( t )  -  H( t-T) J . 


(83) 


It  should  be  noted  that  the  time  tp=(a/c^(r/a  -2sin?/2)  corresponds  to 
arrival  of  the  reflected  pulse  at  the  point  (r,$>).  Let  w  =  w/tu0  then 
for  t^  <  t  <  tj  +  T 


The 


FG(r,t,->)  =  /  pG(r)h(  t-x  )dx 
o 

first  term  of  Eq.  (82)  when  substituted  in  Eq. 


(84) 


(84)  gives: 


P0  /  2sinc?/ 2  t 

I,  - - / - - -  (  c  cosu>Td(2asin4  +  ct-ci-r)dT 

1  4tt  r/a  ‘  2 


1  o 
4tt 


2sin}/ 2 


r/a 


cos[(~)  (2sin4  +  ^^)] 


This  term  would  be  equal  to  zero  for  t  >  tj  +  T.  The  2nd  terra  of  Eq. 
(82)  which  is  of  the  torn  h( r , t , ; )u( t-t { )  when  substituted  on  Eq .  (84) 


leads  to: 


J-2  -  I  coswxh(r,t-x,<{>)  [u(x)-u(x-T)  ]u(t-x-tj)dx 
o 

t-ti  t 

=  /  coscuTh(r  ,t-x  ,rj)dT  =  /  cosw( t-x )h(r ,x ,? )dx 
0  tl 

if  t  <  +  T 

T  t 

=  /  coswTh(r  ,c-t  ,i{i)dx  =  /  cosw(t-x  )h(r  ,i  ,<p  )dx 
o  t-T 

i_f  t  >  tj  +  T 

Substituting  for  h(r,t,<}>)  one  obtains  expressions  for  I2  during  the 
tine  interval  t  <  tj  +  T  as: 


o 

^  4tt 


r/a 


< 2xt i ) (— )  Z  in  e~ir,n(  a  +  2sin2) 
a  n 

n 


C  innx  u  . 
J  e  11  cos(—  x) 

0  X 


dx 

<f> 


PQ  /2sin$/  2  Jln  -  —  — 

=  ~  47  * — 77 —  (2Tfi)  Z  - - -  {inncos(^-  o)  +  ^  sin(~  0) 

“  ,r,  ^  a  a  x 

(-)  -  n 


-inno-, 
-irin  e  n  } 

The  geometrically  reflected  pulse  can  finally  be  expressed  as: 


FG(r,t,o)  =  -  7—  n-~  —  { cos (*—  0)  +  2" i  Z  Zn( - i- 


4tt  '  r/a  A  - 7  Unncos(- 

*  n  * 

(— )  -  n 


+  (“)  sin(—  a)  -  in„  e  inn= 


(85) 
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where  a  =  ct-r/a  +  2sin<}>/2.  It  can  be  observed  that  in  the  case  of 
w/A  =  hn,  the  corresponding  term  in  the  series  summation  would  be 
indeterminate  which  nevertheless  has  a  finite  limit. 


inncT/a  .  ,  -inno  (86) 

e  n  -  mni  e  n 

The  factor  e“inna  represents  an  exponentially  decaying  function  of  time 
in  this  equation. 

It  should  be  observed  that  2n  in  the  above  equations  is  a  function 
of  both  the  radius  of  the  cylinder  (because  of  the  factor  X  =  (c/a)/wc 
in  the  expression  of  the  impedance  function)  and  the  observation  angle 
(y  =  ncos;j>/2).  For  any  values  of  the  dimensionless  parameter  A  and  the 
observation  angle  Eq.  (75)  should  be  solved  numerically  for  the 
complex  value  of  nn  in  the  lower  half-plane.  Eqs.  ( 78 ) — ( 81)  are  then 
used  to  compute  the  corresponding  values  of  £n  which  are  then 
substituted  in  Eqs.  (85)  and  (86)  to  compute  the  geometrically 
reflected  pressure  as  a  function  of  time  and  observation  point 
coordinates.  On  the  other  hand,  changing  the  parameter  to  =  m/w0  in  Eq. 
(85)  will  demonstrate  the  dependence  of  the  scattered  pressure  on  the 
characteristics  of  the  impedance  boundary  as  was  discussed  for  the  case 
of  a  plane  boundary. 
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3.6  Creeping  Waves 

For  an  observer  in  the  farfield,  the  creeping  wave  part  of  the 
scattered  pressure  was  shown  in  Eq.  (53)  to  be: 


One  raay  express  the  Hankel  function  in  terns  of  a  function  related  to 
Airy's  function  as  follows: 


p  »  1 


where  A(q)  is  related  to  the  Airy's  function  and  is  defined  by  [18] 


A(q)  =  -j  /  e*^1  T  ^dx  =  /  cos(qx-x3)dx  , 


A  (q)  +  -9-  A(q)  =  0. 


Differentiating  Eq.  (91)  with  respect  to  the  argument  one  obtains: 

^  H<U(p)  =  H^(1)(p)  ~  -  (|)(|)2/3  e  3  A'  ( q )  (94) 

and  differentiating  Eq.  (91)  with  respect  to  order  one  obtains: 


37  »>>  -  <!><7)2/3  **<0  -  - 


In  addition,  differentiating  Hy  (p)  in  Eq.  (94)  with  respect  to  y  one 
obtains : 


^H.;(1)(p)  ~  (|)(|)  A"(q)  . 


Substituting  Eqs.  (91)  through  (95)  into  Eq .  (90),  expressions  for 
the  functions  Rv(n)  and  Qv(n)  are  obtained  as  follows: 


LIT  '(n)  +  iZ(n)  Hy'  (n)  ] 


2/3  -  jP-i  ( ( 

{(£)(-)  e  3  A'(q)  +  iZ(n)(|)(|)A"(q)} 

Tt  n  TT  n 

Y=nr( 


-1/3  - 

(|)(|)A'(q  ){(|)  e  3  +  iZ(n)(-  f)  ^  } 

n  %  n  A  (q) 

4  Y-Y. 


On  the  other  hand  using  Eq.  (88),  one  can  simplify  this  expression 


further ,  i.e. , 


in  _  2xi 

(I)(i)1/3  e  3  A(q4)  +  iZ(n)(-  |)(|)  e  3  a’  (q^)  =  0 
n 

or  A(q  )  ,  1/3  . 

-X  -  Z(n)(|)  e,l/6 

A  (,,) 


subsituting  for  A(q^)/A  (q^)  in  the  expression  for  dQy  (n )/dY: results 

l 


dQy  (n) 

- -  ~  (~)(~)a' (q5 )  x(n)  , 

dY„  “  "  ^ 

.  2Tti  1/3  ~ 

,  —1/3 - q .  2r.w°\  3 

where  x(^)  =  (~)  e  ~  V- ^ 


,  f6 v  ,  v  -in/3 

and  q „  =  (— )  (Y„-n)e 


Similarly  the  numerator  of  Eq.  (90)  can  be  written  as: 


R  (n)  =  ( h£ 2 ^ ( n )  +  iz(n)  H'(2)(n)j  =H(2)(r,){l  +  iZ(n)  -f 


Y  Y 

Y=Y, 


H*(2)(n) 


Using  the  definition  of  the  Wronskian 


W(H<U(n),  H^2)(n))  =  , 


(100) 


and  combining  Eqs.  (88),  (99)  and  (100)  the  denominator  becomes: 

R  (n)  -  — rr> 

h  "n  H<U(n> 

' i 


Substituting  for  IT  '(n)  in  terras  of  Airy's  function  one  gets: 


R>„(n)  ~ 


AZ(n) 


1/3  -  — 

(1,(1, 1/3  <=  ^  A(,t) 

ti  n 


(101) 


Replacing  R  (n)  and  dQ  (n)/dy.  in  Eq.  (96)  and  (101)  into  Eq.  (90), 
Yl  Yl  1 


R  <n) 
'z 


Kn)  a1"73 

1/3  1 

(6)(6)  A(q£)A  (q£)x(n) 

n 


which  may  be  put  in  final  form  by  substituting  for  A  (q^),  resulting 

in  an  expression  for  K  in  terms  of  .Airy  function: 

Y 

l 


K  (n)  ~  (It—)!  6glTT/2 
*Z  (^)  x(n) 


(102) 
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Substituting  the  asymptotic  approximation  of  Hankel  function, 
valid  in  the  farfield: 


H(l)(A 
Yt  a 


-f- 


r 

IT— 

a 


.  ,  r  it  tt 

l(nI  -'iTi1 


-  »  1 


then  Eq.  (87)  for  the  creeping  wave  part  of  the  scattered  pressure  can 
be  written  as: 


pC(r,t,£) 


Z 

x=  ± 


00  CO 


s=0  i=Q 


/ 


nr/a 


in  .  K  (n) 

4  /  —~— 

-oo  /  n 


if  (n) . 

e  dn  , 


(103) 


where  f(n)  =  y  .(<(>.  +  2tts)  -  fin,  ft  =  C-t-— 

X  A  3 


(104) 


The  integral  in  Eq.  (103)  can  be  computed  by  the  stationary  phase 
technique  but  first  it  is  necessary  to  express  y£  as  a  function  of  n« 
For  that  one  can  use  the  following  approximate  expression  for  Hankel 
function  when  |p|  »  1 


„(1),  x  /  2  e'11T/3  ,  ,1/2  „(1)  rl  ,  3/2 , 

\  (p)  ~  /  j  — JIT  vY-p)  H i / 3  >  . 


(105) 


where  t  =  /3  ,  |-“i  «  1  > 


(n)  =  2i(0/3)  Ai (C  e_TTl/3), 


2/3 


A(q) 


,1/3 


Ai(- 


,1/3 


)  r.) 
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and  Ai(x)  ■  the  Airy  function  [ 1 8 J .  The  asymptotic  expression  for 
H*(D(p)  can  also  be  expressed  as: 


'(1),  .  2  eiTT^3  ,  ,  -iTt/2  u(l)  r  1,  3/2, 

(p) - — - —  (Y-p)  e  H^/3  l-  •3(-2t)  } 

/3 


(106) 


Substituting  Eqs.  (105)  and  (106)  in  Eq.  (88)  results  in  a 
transcendental  equation  on  Y  '• 


— •  3/  2 

3  H<}>  (j  (-  2ts)  } 

Hl/3  {3  (~  2Ts}  } 


6  /-2x, 


(107) 


where  6  =» 


.  1/3 

tp 


(108) 


The  roots  of  Eq.  (107)  can  be  expressed  as  (191 


Tl  3  T£o 


-  5  -  i  T«.  <5 3  +  —  £ 4  -  4  T  c  n  5  5  +  ... 


if  lot  «  1  (109) 


H  =  T; 


JL _  I  _  _ 

2x^00  >5  8x 


1 _  i  + 

-r—  '  •  < 


if  16  I  »  1 


(110) 


,  1.,  ,2/3  in/3 

where:  xia  =  -y(3x£)  e 


x £  are  the  roots  of  J^,j(x) 


J_1/3(x)  =  0 


l._  ,2/3  in/3 

X£oo  =  — (  3x£>  e 


x  £  are  the  roots  of  ^^(x)  ~  J  j  ^(x)  =  0 


The  values  of  the  first  few  roots  of  the  above  equations  are  listed  in 
Table  1. 


Table  1 


t£o  =  uJ£o  e 


iir/3 


=  e 


irr/3 


u0’0 

= 

1.856 

u0 

= 

.808 

U1 .0 

= 

3.245 

U  i  ,00 

= 

2.577 

u2  >  0 

= 

4.382 

u2  »°° 

= 

3.824 

u3  »0 

= 

5.386 

u3 »°° 

= 

4.892 

u4  »0 

= 

6.305 

, co 

= 

5.851 

u6*0 

= 

7.161 

u6  >°° 

= 

6.737 

Another  fora  of  Eq.  (107)  in  terms  of  Airy's  function  can  be  stated 
as : 


-ni/3  '  v 

_e _ A  (q)  =  _  J_ 

A(q)  5 


(111) 


The  roots  of  this  equation  can  be  expressed  in  terras  of  i£  as 
q £  =  6ly/3  ancj  T„  is  given  by  Eq.  (109)  or  Eq.  (110)  depending 

on  the  magnitude  of  the  factor  3  =  Z(p)/ip^3> 
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Substituting  the  impedance  function  Z(p)  into  Eq .  (108)  it  can  be 
observed  that  for  large  value  of  p 

j 6 |  ~  ai  Ip2/3  »  1 

indicating  that  Eq.  (110)  is  the  proper  relation  describing  the  roots 
of  Eq.  (88).  Finally  substituting  for  t£  in  the  right  hand-side  of 
Eq.  (110)  one  gets 


Y-P 

1/3 

P 


T£oo 


_1 _ 

2T£oo 


(112) 


Using  only  the  first  term  of  the  above  series,  a  relation  between  Yjj 
and  n  can  be  written  as 
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It  should  be  noted  that  the  saddle  point  is  a  function  of 

parameters  £,s  and  X.  On  the  other  hand  since  i!  =  (ct-r)/a,  the 

location  of  the  saddle  point  would  also  be  tine  dependent.  For  a  given 

observation  point  (r,^),  the  arrival  time  for  the  creeping  wave  would 

be  t  >  (? ' '  +  2tts)  a/c  +  r/ c  which  implies  that  (i  1  -<j>  ’ '  -2tts)  >  0. 

X 

This  indicates  that  the  saddle  point  falls  on  the  positive  side 
of  the  imaginary  axis  in  complex  n-plane. 

Using  the  conformal  transformation 

w2  =  f(n0)  ~  f(n),  (116) 


The  integration  on  n  in  Eq .  (103)  can  be  approximated  as: 


[*”  _  eif(n)  dn  -  -2  [-Ky1.—  1  ei  If  (n0)-w2  J  rdn 


/T 


/n  n=n(w) 


(3t>  d" 


-  -2  elf(r'oJ  /  F  (n)  e  ^  dw  017) 


where 


K  (n) 


F  (n)  = 


Z  dn 


Expanding  F  (n)  in  Taylor's  series  about  n=no(w=0)>  then 


dF^  (n) 

Fy  (n)  =  F  (n0)  +  [ — - - f  (n  -  n3)  +  .. 

1  l  Hq 


K  (n)  K  (n) 

' Z  /dn N  ,  d  ,  Z  dn ,  , _ 

- (' 

lo  no  /n 


Sn{ 


0> 


(118) 
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From  Eq.  (116)  expanding  f(r\)  in  Taylor's  series  about  no 


w2  =  f(n0)  ~  Lf(n0)  +  f  (n0)(n  -  n0)  +  —  t  (n0)(n  -  n0)  +  •••) 


=  [77  f  (n0)(n  -  n0)2  +  4r  f  (n0)(n  -  n0)3  +  ••• 


(119) 


Differentiating  Eq.  (116)  with  respect  to  w 


n  dw  c  /  \  dr|  2w 

2w—  =-f  (n)  or  —  - ; - 

dn  dw  , '  . 

f  (n) 


(120) 


and  substituting  for  w  in  terms  of  n  we  finally  find 

f  '  '  '  ( n  ) 

Lim(— )  =  Lim - 7 —  [-  -57-  f  (n0)(n-n0)2  -  — -  (n-n0)3  + 

dw  f  (n) 


or  n+nr 


/  f  (n0)  (n-n0) 

-2/ - - -  iim  — 7 - 

z  f  (n) 


-2f  (n0)  (  ~ n — ) 

f  (n0) 


Lim(j£)  =  -  /  — tt — 
dw  £  (n0> 


(121) 


Using  Eqs.  (117),  (118)  and  (121),  the  first  term  of  the  asymptotic 
series  approximation  could  be  obtained  as  follows: 


■h>  KY,(n)  .  . .  .  . ,  n  ‘Sr ,  ^ n  ^ 

;  -  elt(n^dn - 2  elt(n°  - - 

-°®  /  n  0 


> 

‘~tw 

e  dw 


(n0)  o 


_  Kf  ^ 

/  2n  i 

=  /  — n -  - - 

f  (  Do  )  ^  ~  j 


i  f  ( n  , ) 


(  122) 
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Thus  the  creeping  waves  in  Eq.  (103)  are  evaluated  by: 


pc(r,t,<i>)  ~  Z 
A*  ± 


OQ  oo 

z  z 

s=0  £.®0 


e-in/4 


K  (h0) 

»  0 

[— — 


/2iTi 

/  — n - 

f  <n0) 


eif(n0) 


] 


(123) 


The  terms  if(n0)  and  /2ffi/f ' ' (n0)  which  are  needed  in  Eq.  (123)  can  be 
obtained  by  using  Eqs.  (114)  and  (115) 


if(n0) 


=  -  2( 


aiia> 


3/2 


i«  3/0 

+  2ns  r/z 

in  -  (4»x  +  2ns)J1/Z 


and 


/ 


2ni 
•  » 


f  (n0) 


3/n 


• »  3/4 

lu^<0»(<Jix  +  2tts)) 

_  FT  a 

[3(n  -  -  2its)r/4 


(124) 


Finally  Eq.  (123),  representing  the  creeping  wave  contribution  to  pulse 
scattering  can  be  rewritten  as 

r—  OO 


pC(r,t,$) 


/  —  Z 
7T—  A*  ± 


z 

s=0 


Px,(r.e,«*(etx,) 


where 

PAs  <*•*•♦> 


_ _  K  CnQ) 

t  -li _ 

&  f  (Hq)  / n  q 


e  if(n0) 


(126) 


( 127-a) 


The  unit  step  function  H(t-txs)  represents  the  arrival  of  the  creeping 
wave  corresponding  to  known  values  of  A  and  s 

-  (}”+  2.s)(f)  (127-b) 

The  pulse  response  as  expressed  by  Eq.  (126)  can  be  substituted  in 
Duhamel's  integral  to  find  the  scattered  pressure  due  to  creeping  waves 


only  for  any  incident  plane  wave.  Letting  the  time  signature  of  the 
incident  pulse  be  described  by: 


Finc(t)  -  cos(ut) (H(t)  -  H(t-T) J  (128) 

then  for  tXs  <  t  <  t  <  tXs  +  T,  the  creeping  wave  pulses  are  given  by: 

c  /2 —  “  t-txs 

Fc(r,t,4>)  *  /  —  Z  Z  /  cosCojtJP^  (r ,t-x  ,<j>)dT 

r  <  l  «  as 

IT—  A**  i  S=0  0 

a 


“  t 

z  z  /  cos[w(t-T)]p  (r  ,T  ,<J>)dT 

X-  ±  s=0  t,  Ks 

Xs 

(129) 


After  a  change  of  variable  x  =*  ct/a  and  substituting  for  u>  in  terns  of 
dimensionless  parameters  ui  =  (o>/X)(c/a)  one  finds: 


— .  ct/a  —  , 

FC  (  r ,  t ,  $  )  =  /  —  Z  Z  /  cos((“)(£E-x)jp  (r,x,*)£* 

tt!  X  s  ctXs/a  X  a  AS  c/a 

Finally  multiplying  by  the  factor  (-P0C/4ir),  which  was  dropped  out 
earlier  one  obtains  the  time  signature  of  the  creeping  wave  pulse  as: 

FC(r,t,<?)  *  -  /% 

TT— 

a 

in  terms  of  the  new  variable  x,  2  =  x  -r/a.  The  integral  in  Eq.  (130) 
can  be  computer  numerically  for  any  non-dimensional  time  ct/a. 

For  t  >  tXs  +  T 

FC(r,t,$)  «  -  /  ~ 

TT— 

a 

(131) 


4tt 


ct 

CO  -  - 

„  r  a  f/u>.,ct 

Z  J  cos[(-)( — 

s-0  c(t-T)  X  a 

a 


x)}pXsU,xA)d* 


ct/a  — 

-T-  L  Z  J  cos[(^)(^~  -  x)]pXs(r,x,(|))dx 

X-±  s-0  ctXs/a  X  (130) 
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Ic  should  be  noted  that  in  computing  p  (r , t , <J> )  according  to  Eq.  (127), 

As 

a  summation  on  index  Z  should  be  carried  out,  but  the  factor 
represents  an  exponential  decay  with  Z  (see  Eq.  (124),  so  considering 


CHAPTER  IV 


NUMERICAL  RESULTS  AND  CONCLUSIONS 


4.1  Introduction 

Numerical  evaluations  of  the  scattered  field  for  various  types  of 
incident  plane  pulses  are  presented  in  this  chapter.  Creeping  waves 
and  geometrically  reflected  waves  are  plotted  separately.  The 
parameters  are  varied  to  explore  the  dependence  of  the  scattered 
pressure  on  the  observation  angle  9  and  impedance  covering  Z.  All 
calculations  were  made  for  an  observer  in  the  far  field.  The  computed 
time  signatures  describe  the  ratio  of  the  amplitude  of  the  scattered 
pressure  to  that  of  the  incident  pressure,  plotted  versus  non- 
dimensional  time  ct/a. 

4.2  Impedance  Model 

Since  the  scattered  field  is  a  function  of  the  impedance  coverage 
as  well  as  the  radius  of  the  cylinder,  one  can  plot  the  pressure  at 
specific  values  of  parameters  aj,  cj  and  A.  Another  factor  which 
influences  the  total  time  signature  is  the  parameter  w.  Once  the 
impedance  function  and  observation  angle  are  fixed,  Eq.  (75)  is  solved 
numerically  for  the  poles  nn  in  the  i°wer  half  of  the  complex  plane. 
Tables  2  and  3  list  some  of  these  poles  corresponding  to  A  *  .5  and 
A  *  2.0  , respectively .  The  numerical  solution  of  Eq.  (75)  has  been  done 
with  the  help  of  1MSL  zero  finder  subroutine  which  works  for  an 
analytic  function  of  complex  variables. 

4.3  Geometrically  Reflected  Pressure 


The  first  three  roots  of 


H^(n)  +  iZ(n)  H  ^(n)  =  0,  y  ■  ncosd/2 

with 


*° 

X  =  2.0, 

nl 

aj  =  2.0  and  bj  =  1.1 

n2 

n3 

0 

-.2340,-. 1745 

10 

-.2342,-. 1748 

20 

-.2348,-. 1763 

30 

-.2357, -.1790 

-1.3990,- 

.0000 

40 

-.2371, -.1827 

-.5350, -.0478 

-1.4300,- 

.0000 

50 

-.2387,-. 1875 

-.5400, -.0001 

-1.4900,- 

.0001 

60 

-.2406, -.1903 

-.5501, -.0300 

-1.5410,- 

.0005 

70 

-.2425, -.2004 

-.5741, -.0324 

-1.6180,- 

.0094 

80 

-.2445, -.2076 

-.5957, -.0270 

-1.7150,- 

.0260 

90 

-.2464, -.2158 

-.6233, -.0226 

-1.8350,- 

•  0b07 

100 

-.2483, -.2244 

-.6507, -.0209 

-1.9753,- 

.1314 

110 

-.2561, -.2333 

-.7039, -.0251 

-2.0918,- 

.2806 

120 

-.2518, -.2423 

-.7607, -.0414 

-2.0560,- 

.4126 

130 

-.2532, -.2510 

-.8264, -.0826 

-2.001  ,- 

.3662 

140 

-.2544, -.2590 

-.8779,-. 1623 

-2.0540,- 

.1817 

150 

-.2554, -.2661 

-.8804, -.2525 

160 

-.2561, -.2717 

-.8554, -.3111 

170 

-.2565, -.2754 

-.8331, -.3407 

-3.7821,- 

.0001 

180 

-.2567, -.2770 

-.8231, -.3512 

-3.8990,- 

.3410 

The  oscillatory  nature  of  the  GRP  for  an  incident  plane  wave  pulse  is 
due  to  the  interference  of  the  geometrically  reflected  my  and  the 
diffracted  rays  fron  the  shadow  boundaries.  A  second  type  of  an 
incident  plane  wave  considered  here  is  a  tine  harmonic  pulse  train. 
Figures  25-29  represent  the  three  dimensional  plot  of  the  scattered  GRP 
versus  time  and  observation  angle  varied  from  0°  to  180°.  It  can  be 
seen  that  the  peaks  at  each  observation  angle  correspond  to  a  specific 
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arrival  time,  which  is  related  to  that  observation  angle.  The 
amplitude  of  the  scattered  GRP  on  the  other  hand  shows  a  distinct 
dependency  on  observation  angle.  Figure  30  demonstrates  a  two 
dimensional  plot  of  the  back  scattered  GRP  for  the  same  incident  plane 
wave  pulse  train.  Since  the  duration  of  the  incident  pulse  is  limited 
to  0  <  u)t  <  Sir ,  the  arrival  time  for  the  scattered  field  is  described 
by: 

(r/a  -  2sin$/2)  <  ct/a  <  8n(X/w)  . 

The  GRP  after  the  termination  of  the  pulse  would  be  an  exponentially 
decaying  function  as  demonstrated  by  negligible  amplitudes  in  Figure 
30.  Figures  31  and  32  exhibit  the  significance  of  the  nondimensional 
frequency  of  the  incident  pulse  (u  on  the  back  scattered  GRP.  It  is 
observed  that  unlike  an  infinite  plane  boundary,  there  are  little 
change  in  the  amplitude  of  the  reflected  pressure  vs.w.  Comparison 
of  Figures  30  and  33  shows  the  backscattered  GRP  for  cylinders  with 
different  nondimensional  radii  X.  The  change  in  X  from  Figure  30  to 
Figure  33  by  a  factor  of  four  corresponds  to  reducing  the  radius  by 
1/4.0.  Since  the  ratio  of  r/a  is  kept  constant,  any  comparison  between 
Figures  30  and  33  is  meaningful  only  if  the  observation  point  is  kept 
fixed  in  space.  This  indicates  that  the  time  signature  as  shown  in 
Figure  33  should  be  multiplied  by  the  factor 

/10/(9X2/Xl+1)  =  .51 

so  that  the  geometrically  spreading  factor  is  eliminated. 
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Figures  34  and  35  represent  the  backscattered  field  for  FM  plane 
pulse  incidence.  Once  again  these  figures  do  not  show  a  significant 
variation  in  the  amplitude  of  the  backscattered  GRP  with  changing 
frequency,  a  finding  which  is  in  agreement  with  those  of  Figures  30  and 
32. 

Finally  Figures  36  and  37  show  polar  plots  for  the  peak  amplitude 
of  the  GRP  time  signature  as  a  function  of  the  observation  angle  for  a 
time  harmonic  pulse  train.  The  forward  scattering  is  zero  as  was 
stated  in  Eq.  (86). 

4.4  Creeping  Waves 

We  saw  that  the  fronts  of  the  creeping  pulses  move  with  the 
free-space  sound  velocity  and  propagate  around  the  shadow  side  of  the 
cylinder  back  to  insonified  region,  and  then  circle  the  cylinder  an 
infinite  number  of  times.  The  exponential  factor  eif^r|o^  in  E9S*  024) 
and  (126)  indicates  the  fast  decaying  nature  of  these  pulses  as  they 
circumnavigate  the  cylinder's  surface.  Figure  38,  39  and  40  show  the 
pulse  response  for  the  first  mode  of  the  circumferential  wave  (s=0)  at 
<J>  *>  180°,  90°  and  0°  respectively.  The  discontinuity  in  Figure  39 
corresponds  to  the  two  different  arrivals  of  the  creeping  waves  from 
the  two  different  directions.  It  should  be  pointed  out  that  the  second 
modes  of  the  creeping  waves  (s=l)  are  so  highly  attenuated  that  they 
show  no  visible  discontinuities,  unless  plotted  separately  with  a 
different  scale  factor.  Even  for  the  first  creeping  wave  mode,  as  they 
travel  from  0°  to  180°,  the  pulses  seem  to  attenuate  significantly. 

The  first  mode  of  the  creeping  waves  for  a  tine  harmonic  pulse  train  is 
piotted  in  Figures  41-44  for  different  observation  angles.  All  these 
figures  indicate  that  the  amplitude  of  the  creeping  wave  increases  as 
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the  observation  angle  changes  from  180°  to  0°,  i.e.,  most  of  the 
creeping  wave  field  is  concentrated  in  the  shadow  zone  of  the  cylinder. 
Figures  45  and  46  correspond  to  the  second  arrival  of  these  waves  at 
180°  and  0°,  respectively.  The  first  arrival  for  the  FM  pulse  train  is 
shown  in  Figures  47  and  48.  The  peak  amplitude  of  the  creeping  waves 
time  signature  is  plotted  in  polar  coordinates  in  Figures  49  and  50. 
Since  these  figures  are  plotted  for  r/a  a  10  but  for  different  values 
of  the  nondimens ional  radius  of  X,  the  amplitude  distribution  as 
presented  in  Figure  50  should  be  multiplied  by  the  factor 

✓l0/(9l2/X1  +  1)  -  0.51 
to  facilitate  comparison  in  Figures  49  and  50. 

4.5  Conclusions 

The  prediction  of  the  impulse  response  of  an  impedance  covered 
cylinder  was  the  main  objective  of  this  study.  The  scattered  pressure 
was  obtained  for  an  incident  plane  wave  pulse,  a  time  harmonic  pulse 
train  and  FM  pulses.  The  computations  were  carried  out  for  an  observer 
in  the  far  field  and  at  different  observation  angles.  The 
nondimensional  frequency  and  radius  w  and  A  were  varied  in  order  to 
assess  the  influence  of  the  frequency  of  the  incident  wave  and  the 
radius  of  the  cylinder  on  the  scattered  pressure.  The  following 
conclusions  can  be  drawn  from  the  study  of  the  analytic  solutions  and 
their  numerical  evaluation: 

a.  The  directivity  of  the  peaks  of  the  geometrically  reflected 

pressure  time  signature  shifts  towards  the  shadow  zone  as  A  is 
increased  (or  as  the  radius  is  decreased)  with  new  diffraction 
peaks  appearing  and  moving  towards  the  shadow  zone. 
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b.  In  contrast  to  the  infinite  flat  plane  boundary,  changes  in 
the  nondimensional  frequency  w  do  not  affect  the  backscattered 
field  significantly  but  it  changes  the  directivity  of  the 
peaks  of  the  scattered  field. 

c.  The  directivity  of  each  peak  of  the  creeping  waves  is  close  to 
being  unidirectional. 

In  the  numerical  evaluations,  the  choice  of  only  s  =  0  and  l  * 
4  would  incorporate  the  major  contribution  arising  from  the 
infinite  series  summation. 

d.  The  forward  scattering  of  the  circumferential  creeping  waves 
could  be  increased  by  decreasing  the  nondimensional  radius  X 
or  increasing  a.  This  agrees  with  the  preceding  studies  for 
the  scattering  from  a  rigid  cylinder. 

4.6  Suggested  Future  Research 

The  analytical  solution  of  the  problem  was  based  on  the  separation 
of  variables  and  the  Sommerf ield-Watson  transformation.  This  technique 
could  be  applied  to  other  simple  geometries  such  as  spheres  and 
spheroids.  The  Watson  transformation  allows  decoupling  of  the 
geometrically  reflected  field  from  the  creeping  waves  which  is 
essential  in  the  understanding  of  the  multiple  echoes  that  return  in 
the  case  of  scattering  of  a  simple  plane  wave  pulse.  Other  techniques 
such  as  the  T- Matrix  method  or  finite  element  method  could  be  applied 
to  solve  for  the  total  scattered  field.  The  impedance  model,  further, 
can  be  extended  to  multi-degree  of  freedom  oscillators,  thus  creating 
new  resonance  frequencies  for  the  impedance  boundary. 
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